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20.  ABSTRACT  / 

This  document  is  the  final  report  for  the  Office  of  Naval  Research/ 
grant  number  N00014-79-C-0624,  entitled  "Parametric  Acoustic  Array  Formation 
via  Weak  Collinear  and  Noncollinear  Interaction  in  Dispersive  Fluids."  -*This 
final  report  contains  a  historical  review  of  the  technical  and  administrative 
highlights  of  the  research  project. The  details  of  the  scientific  findings 
are  available  in  a  separate  technical^ memor andugof  the  same  title,  available 
through  the  Defense  Technical  Information  Center/^  Briefly,  the  project  was  a 
theoretical  investigation  of  the  effects  of  dispersion  on  the  nonlinear 
formation  of  a  difference-frequency  beam  from  two  primary  beams.  The  approach 
involved  both  analytical  and  numerical  solutions  to  the  nonlinear  and  dis¬ 
persive  paraxial  wave  equation  with  emphasis  on  Gaussian  beam  solutions.  In 
the  case  of  collinear  primary  beams,  known  results  were  reproduced,  but  with 
some  significant  errors  corrected.  Noncollinear  interaction  of  finite  width 
beams  was  investigated  for  the  first  time  in  this  project.  Some  unexpected 
results  were  discovered  in  this  case:  it  was  found  that  the  finite  beam  width 
caused  the  difference  beam  to  have  a  directivity  maximum  near  the  axis  of  the 
primaries,  rather  than  in  the  direction  of  the  primary  wave-vector  difference, 
and  decreased  the  conversion  efficiency  to  significantly  below  the  plane-wave 
efficiency.  Further  research  in  this  area  is  indicated. 


neaTiew  e*  tw m  e 


ONR  PINAL  REPORT 
N00014-79-C-0624 

1  August  1979  -  30  September  1983 


Parametric  Acoustic  Array  Formation  via  Weak  Collinear  and  Noncollinear 
Interaction  in  Dispersive  Fluids:  Final  Report 


Julian  D.  Maynard 
The  Pennsylvania  State  University 
University  Park,  PA  16802 


Reproduction  in  whole  or  in  part  is  permitted  for  any  purpose  of  the 

United  States  Government 


This  final  report  presents  a  historical  review  of  the  technical  and  admin¬ 
istrative  highlights  of  the  research  program  undertaken  under  ONR  Grant 
N00014-79-C-0624.  Details  of  the  scientific  findings  can  be  found  in  the  technical 
report  "Parametric  Acoustic  Array  Formation  via  weak  Collinear  and  Noncollinear 
Interaction  in  Dispersive  Fluids,"  written  by  Mark  F.  Hamilton  who  was  the 
graduate  student  supported  by  the  grant*  This  technical  report  may  be  referenced 
through  the  Defense  Technical  Information  Center  Report  No.  ADA  130533,  or  through 
the  Applied  Research  Laboratories,  University  of  Texas,  Austin,  Report  No. 
ARL-TR-83-19. 

The  original  idea  for  this  research  was  formulated  by  Francis  H.  Fenlon, 
Associate  Professor  of  Engineering  and  Senior  Research  Associate  of  the  Applied 
Research  Laboratory  at  The  Pennsylvania  State  University.  Briefly,  the  idea  was 
to  study  the  effects  of  induced  and  inherent  dispersion  on  the  conversion-efficiency 
of  an  acoustic  parametric  array.  In  a  parametric  array,  two  primary  sound  fields 
with  frequencies  u>i  and  u> 2  and  wave  vectors  kj  and  k2  interact  in  a  nonlinear 
medium  to  produce  secondary  fields  at  sum  and  difference  frequencies.  If  the 
medium  is  nondispersive  (the  speed  of  sound  c  is  independent  of  frequency)  then 
the  nonlinear  conversion  is  most  efficient  for  collinear  plane  waves,  since  in 
this  case  the  phase  of  the  secondary  wave  matches  the  beating  of  the  primary  waves 
[i.e.,  k^  =  (go ^  ±  GJ2)/c  »  k }  t  k2 1  •  However  if  the  medium  is  dispersive,  then 
maximum  efficiency  occurs  when  the  primary  and  secondary  waves  intersect  at  a 
particular  angle,  and  only  one  type  of  secondary  wave  [sum  or  difference 
frequency]  is  efficiently  generated  depending  on  whether  the  dispersion  is 
positive  [c  increases  with  frequency]  or  negative  [c  decreases  with  frequency]. 

In  sonar  applications  the  difference  frequency  wave  is  of  greater  importance, 
and  the  use  of  dispersion  to  enhance  its  generation  is  of  practical  interest. 
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Also  of  practical  interest  is  the  effect  of  the  finite  size  of  the  projector 
which  generates  the  primary  beams;  the  finite  projector  size  gives  the  primary 
waves  an  extended  angular  spectrum  (referred  to  as  beam-spreading  or  diffraction) 
so  that  their  interaction  is  intrinsically  noncollinear  and  so  that  the  finite 
wave  amplitude  decreases  from  geometrical  spreading  as  well  as  from  thermo- 
viscous  absorption  in  the  fluid.  The  original  proposal  considered  the  effects 
on  seconary  beam  generation  from  two  sources  of  dispersion:  artifical,  boundary- 
induced  dispersion,  and  inherent  dispersion  arising  from  relaxation  mechanisms. 

During  the  course  of  the  research,  emphasis  shifted  solely  toward  the  latter 
source  of  dispersion.  Although  there  exist  papers  dealing  with  the  effect  of 
dispersion  on  infinite  plane  waves,  the  results  in  Mr.  Hamilton's  report  are 
the  first  to  describe  dispersion  effects  with  finite  beams.  His  results 
uncover  some  unexpected  behavior  and  indicate  that  further  research  should  be 
undertaken;  his  work  also  reveals  many  errors  in  the  existing  literature. 

These  results  will  be  discussed  briefly  in  this  report. 

Mr.  Hamilton  began  working  on  this  project  under  the  guidance  of  Professor 
Fenlon  in  1980,  and  after  a  year  he  mastered  the  relevant  concepts  existing  in 
the  literature  and  began  working  on  new  calculations.  Initial  results  were 
reported  at  the  Ottawa,  Canada  meeting  of  the  Acoustical  Society  of  America  [J. 

Acoust •  Soc.  Am.  69^,  S81  (1981)]  in  May  1981.  In  June  1981  the  theoretical  acoustics 
community  suffered  considerable  loss  when  Professor  Fenlon  passed  away.  Through 
his  diligence  and  excellence  as  a  teacher,  he  had  advanced  the  project  to  a  point 
where  Mr.  Hamilton  could  see  it  to  completion.  Mr.  Hamilton  had  become  suffi¬ 
ciently  versed  in  the  field  that  he  was  able  to  substitute  for  Professor 
Fenlon  in  delivering  the  basic  nonlinear  acoustics  lectures  at  the  Penn  State 
Underwater  Sound  School  in  October  1981  and  again  in  1982. 
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In  order  to  allow  the  research  project  and  its  ONR  funding  to  continuer 
Physics  Professor  Julian  D.  Maynard,  a  friend  and  colleague  of  Professor  Fenlon, 
assumed  the  responsibilities  of  principal  investigator  and  advisor  for  Mr. 
Hamilton's  academic  affairs.  The  transfer  of  the  grant  to  Professor  Maynard  was 
officially  approved  in  November,  1981  as  a  modification  A00002  of  the  original 
proposal.  For  the  advanced  scientific  guidance  of  the  nonlinear  acoustics 
research,  aid  was  most  generously  offered  by  Professor  David  J.  Blackstock  of  the 
University  of  Texas.  In  addition  to  being  distinguished  in  the  fields  of  physical 
and  nonlinear  acoustics.  Professor  Blackstock  was  a  personal  friend  of  Professor 
Fenlon.  In  order  for  the  scientific  advising  to  be  administered  conveniently, 
it  was  necessary  that  Mr.  Hamilton  complete  his  theses  work  at  the  University 
of  Texas;  his  degree  would  still  be  conferred  through  the  Acoustic  Program 
at  Penn  State.  To  qualify  as  an  off-campus  degree  candidate,  Mr.  Hamilton  was 
required  to  complete  his  Ph.D.  candidacy  exam.  In  November  1981  he  passed  this 
exam  with  an  above-average  score,  after  which  his  Ph.D.  degree  committee,  with 
Professor  Blackstock  established  as  an  external  member,  was  formed.  After 
completing  these  academic  requirements,  Mr.  Hamilton  moved  to  Texas.  Reimbursement 
for  this  move,  costing  approximately  $1000,  was  approved  and  paid  through  the  ONR 
grant. 

Because  of  the  delays  in  the  progress  of  the  research  caused  by  the  change  in 
theses  advisor,  the  move  to  Texas,  etc.,  it  was  agreed  that  a  non-cost  extension  of 
the  grant  would  be  necessary  to  complete  the  project.  A  new  termination  date  of 
May  31,  1983  (giving  Mr.  Hamilton  an  additional  ten  months  to  compete  his  Ph.D. 
thesis)  was  requested  and  approved  in  March  1982.  A  revised  budget  including  Mr. 
Hamilton's  travel  expenses  to  the  University  of  Texas  and  his  various  research 
expenses  there,  was  also  approved  at  this  time. 
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In  order  to  conveniently  monitor  Mr.  Hamilton's  progress,  it  was  agree  that 
the  bi-annual  conventions  of  the  Acoustical  Society  of  America  would  be  used 
as  meeting  places  for  Mr.  Hamilton  and  Professor  Blackstock  from  Texas  and  Professor 
Maynard  from  Penn  State.  Such  meetings  took  place  at  the  April  1982  convention  in 
Chicago  and  the  November  1982  convention  in  Orlando,  Florida.  At  the  Orlando 
meeting,  Mr.  Hamilton  delivered  a  talk  [J.  Acoust.  Soc.  Am.  12/  S41  (1982)]  in 
which  he  presented  three-dimensional  plots  of  the  sound  field  resulting  from 
two  collinear,  interacting  Gaussian  beams.  These  clearly  showed  that  a  sufficient 
amount  of  dispersion  caused  the  maximum  radiation  to  be  shifted  off  axis,  at 
an  angle  9  given  by  k^  -  kj  M  k-  cos  0.  An  important  aspect  of  this  work  was 
the  correction  of  numerous  errors  existing  in  the  literature. 

While  at  the  University  of  Texas,  Mr.  Hamilton  completed  the  analytical 
and  numerical  investigation  of  the  collinear  interaction,  and  then  extended 
the  calculations  to  the  case  of  noncolliner  interaction.  The  culmination  of 
this  work  was  Mr.  Hamilton's  thesis  (which  is  reproduced  in  the  technical 
memorandum  mentioned  in  the  introduction  of  this  report)  a  shortened  version 
of  the  thesis  submitted  for  publication  to  the  Journal  of  the  Acoustical 
Society  of  America  (pending),  and  a  paper  prsented  at  the  Eleventh  International 
Congress  on  Acoustics  in  Paris.  An  excellent  summary  of  the  scientific  accomplish¬ 
ments  of  this  project,  including  principal  techniques  and  results,  is  contained 
in  Appendix  A,  which  is  a  reproduction  of  the  paper  prepared  by  Mr.  Hamilton 
for  the  Paris  meeting.  A  significant  result  was  the  effect  of  the  finite  beam 
size  in  the  noncollinear  beam  interaction.  When  the  beams  intersect  at  the 
resonance  angle,  the  position  of  maximum  radiation  was  not  in  the  direction  of 
kj  -  k2#  but  remained  close  to  the  primary  beam  axes.  This  was  a  result  of 
the  directionality  of  the  finite  width  interaction  region.  Another  interesting 


result  was  the  United  secondary  field  generation  caused  by  the  reduced  inter¬ 
action  length  of  the  intersecting,  finite  width  primaries.  Further  investigation 
of  this  compromise  in  the  noncolliner,  finite  beam  situation  is  indicated. 

In  order  to  complete  the  technical  memorandum  and  final  report  covering  the 
project  and  to  include  the  Paris  conference  as  part  of  the  project,  it  was 
necessary  to  request  a  second  no-cost  extension  of  the  grant.  A  termination  date 
of  September  30,  1983  was  approved  in  May  1983. 

On  April  8,  1983,  Mr.  Hamilton  passed  his  comprehensive  examination  with 
above  average  scores,  and  on  the  following  day  he  successfully  defended  his 
thesis,  thus  completing  the  requirements  for  the  Ph.D.  degree. 
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Introduction.  The  effect  of  dispersion  on  parametric  acoustic  arrays 
formed  by  collinear  interaction  of  weak  finite  amplitude  waves  is  well 
understood.1  Inherent  dispersion  in  fluids  results,  for  example,  from 
relaxation  phenomena,  or  inhomogeneities  such  as  bubbles.  Compensation  for 
the  detrimental  effects  of  dispersion  on  the  nonlinear  interaction  of  two 
finite  amplitude  plane  waves  is  achieved  when  the  two  waves  intersect  at  a 
suitable  angle."  Such  an  approach  was  applied  by  Kozyaev  and  Naugol’nvkh3 
in  their  analysis  of  parametric  arrays  formed  by  collimated  plane  primary 
waves  in  bubbly  water.  In  this  paper  we  extend  our  previous  analytical 
work, 1  where  we  used  Gaussian  primary  beams  to  account  for  dissipation  and 
diffraction  as  well  as  dispersion,  to  the  case  of  noncollinear  interaction. 

Theory.  We  begin  by  modifying  Kuznetsov’s4  nonlinear  paraxial  wave  equa¬ 
tion  to  account  for  dispersion.  In  the  frequency  domain  our  modified 
equation  assumes  the  parabolic  form 
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where  p^(x,y,z)  is  the  acoustic  pressure  spectral  amplitude  of  a  quasi¬ 
plane  wave  traveling  in  the  positive  z  direction,  p ^  *  p^  represents  convo¬ 
lution  in  the  frequency  domain,  P0  is  the  ambient  density  of  the  fluid,  cQ 
is  a  reference  sound  speed,  and  3  is  the  coefficient  of  nonlinearity  (e.g., 
S*(Y+l)/2  for  gases,  where  y  is  the  specific  heat  ratio).  Dispersion  is 
taken  into  account  via  the  imaginary  part  of  +  i(kw  -  uj/c0)  ,  where 

is  the  attenuation  coefficient,  and  k^-co/c^  is  the  wavenumber,  where  cw  is 
the  sound  speed  at  angular  frequency  a).  See  Ref.  5  for  a  similar  treatment 
of  the  Burgers  equation. 

By  restricting  our  analysis  to  weakly  nonlinear  interactions,  we  can 
solve  Eq.  (1)  via  successive  approximations  (see  Ref.  6,  Appendix  A). 

First  we  seek  the  solution  for  a  primary  beam  whose  direction  of  propaga¬ 
tion  forms  an  angle  <f>  with  the  z  axis  in  the  x-z  plane.  For  a  Gaussian 
beam  the  boundary  condition  becomes 


Pj(x,y,0)  •  PGj  exp{-(x2  +  y2)/e2  -  iw^x/c^ 
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where  j  ■  1,2  represents  and  PG  is  the  peak  source  amplitude,  and  eQ 
determines  the  spot  size  of  the  beam.  In  the  first  approximation,  where  we 
neglect  p,^  *  p^,  the  solution  for  the  primary  beams  is  thus  found  from 
Eq.  (1)  to  be 
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iw0e£/2c0  is  the  mean  collimacion  distance  of  the 


where  we  have  introduced  the  dimensionless  coordinates  X 
and  Z*z/z0,  where 

primaries,  and  (J.  ■  (u)^  +  0^)  /2  is  the  mean  primary  frequency.  We  also  let 
L-u>j/u>0,  aj*ajz0,  and  5  *  oo0c0/2c0.  In  the  farfield  paraxial  region 
(Z  >>  1)  we  can  let  X  *  DZ9X  and  Y  *  DZ9y*  where  9X  and  9y  are  azimuthal 
angles  in  the  x-z  and  y-z  planes,  respectively,  whereby  the  directivity 
function  for  the  primaries  is  found  from  Eq.  (3)  to  be 
exp{-fi^D2[(9x  -  4>j) 2  +  92]}  . 

To  obtain  a  first  approximation  of  the  solution  for  the  difference 
frequency  signal  (w_  -  -  0J2) ,  we  solve  Eq.  (1)  after  replacing  pfJL>  *  p^  by 

p^pS,  where  p£  is  the  complex  conjugate  of  P2*  Without  loss  of  generality 
we  let  #2  “  0  and  which  creates  the  geometry  depicted  in  Fig.  1. 

Solving  Eq.  (1)  for  p_ 
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where  A  »  1  “  i2Z/Q_)  ,  B  »  Z  +  i(Q_  +  2S2^£22/^— )  *  ft_*w_/u)0, 

PQ  »  8pG^pG2/o,jc|,  and  aj  *  aT  -  12<SJ)_d2,  where  a^  *  ^al  +ot2  “a_)2o*  The  dis¬ 
persion  coefficient  is  <5*1-  |l£,  -  1c2| /k_,  which  for  collinear  interaction 
(i.e.,  0*0)  reduces  to  <$0  *  1  -  (k.  -  kO 7k_,  where  oQ  is  a  measure  of  the 
inherent  dispersivity  of  che  fluid,  we  obtain  $0 *  0  for  no  dispersion, 

60  >  0  for  normal  dispersion  (c^>C2>c_),  and  50  <  0  for  anomalous  disper¬ 
sion  (cj  <  C2  <  O  . 

Results.  We  first  use  Eq.  (4)  to  calculate  the  difference 
frequency  field  for  collinear  interaction  with  S0 *  0,  -0.01, 
and  0.01.  The  field  plots  are  presented  in  Fig.  2,  where  we 
have  chosen  ft_*0.1  (note  that  (2^*l  +  n_/2  and  ^2  ™  1  ~  f2_/2) , 
a>^  ■  0.1,  a_  *  0,  and  D*30.  Since  the  field  is  axisyumetric 
for  0*0,  we  let  ^"X^  +  Y^.  when  there  is  no  dispersion, 
as  depicted  in  Fig.  2(a),  the  virtual  sources  interact  syn¬ 
chronously  with  the  difference  frequency  signal,  giving  rise 
to  the  end-fire  radiation  pattern.  As  seen  from  Figs.  2(b) 
and  2(c),  dispersion  causes  asynchronous  Interaction,  with 
the  interaction  region  being  clearly  defined  by  the  resulting 
spatial  ost  ..latlonr  For  normal  dispersion  (<So>^)i  c^« 
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(a)  6o«0  (b)  5q  *  -0.01  (c)  0Q  *  0 .0 1 

FIGURE  2  -  Pressure  fields  where  $*Q,  ft_  =  0.1,  aT*0.1,  a_  =  0,  D  *  30 


phasing  of  the  virtual  sources  corresponds  to  a  speed  which  is  greater 
than  c_,  while  for  anomalous  dispersion  (6o<0),  the  opposite  is  true. 
Radiation  from  the  interaction  region  is  thus  analogous  to  that  from 
bending  waves  on  an  infinite  plate.  When  5o>0,  the  direction  of  maximum 
radiation  is  shifted  off  axis  to  an  angle  which  approaches 
cos^tCk^  -  ’<2)  /k_]  for  collimated  plane  wave  interaction.  When  5o<0,  the 
radiation  is  of  an  evanescent  nature  and  the  maximum  remains  on  axis. 


Only  for  50 > 0  can  phase  matching  in  the  interaction  region  be 
attained  via  intersection  of  the  primaries.  Phase  matching  is  achieved  by 
increasing  the  angle  of  intersection,  j>9  until  the  wavenumber  of  the  vir¬ 
tual  source  distribution,  |kj-k2L  equals  k_ .  We  will  investigate  the 
effect  of  noncollinear  interaction  by  way  of  the  farfield  radiation  pattern 
of  the  difference  frequency  signal.  The  farfield  asymptotic  form  of 
Eq .  (4)  is  extremely  complicated  and  will  not  be  presented  here.  However, 
tremendous  simplification  results  for  an  array  where  attenuation  restricts 
nonlinear  interaction  to  the  nearfield  of  the  primaries.  Thus,  for  large 
a^.  and  small  the  farfield  asymptotic  form  of  Eq .  (4)  reduces  to 


H2D2P  *a-Z  exo(-(r.2D2/2)[(9  -  (8,  /8  )$) 2  +  32] }  -ifl  D2Z02 
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where  P_ * P_( -x* 3y *z)  and  + 

When  (p  *  0  and  5o*0,  the  denominator 
of  the  directivity  function  in  Eq .  (5) 
reduces  to  Westervelt 's7  result  for 
an  absorption-limited  array. 

To  calculate  farfield  radiation 
patterns  for  various  values  of  we 
evaluate  Eq.  (4)  with  Z*00,  using  the 
same  parameters  as  in  Fig.  2(c).  In 
addition,  it  is  assumed  that  both 
primaries  propagate  at  the  same  speed. 
The  angle  <pQ  for  which  5*0  thus  be¬ 
comes  0.82°.  In  Fig.  3  are  presented 
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(a)  $  *  0°  (b)  5*0.7°  (c)  $*0.82° 


FIGURE  4  -  Farfield  radiation  patterns 

the  farfield  radiation  patterns  for  the  x-z  plane  (i.e.,  where  6V  *  0) . 
Likewise,  complete  radiation  patterns  are  shown  in  Fig.  4. 

From  Figs.  2(c),  3,  and  4(a),  we  find  that  collinear  interaction  with 
5Q >  0  produces  a  dip  in  the  center  of  the  radiation  pattern.  As  0  ap¬ 
proaches  4>o»  the  disappears,  and  the  field  resembles  that  depicted  in 
Fig.  2(a).  Increasing  0  beyond  0O  causes  geometric  dispersion  to  overcom¬ 
pensate  for  the  inherent  dispersivity  of  the  fluid.  The  net  result  of  such 
overcompensation  is  a  field  resembling  that  depicted  in  Fig.  2(b),  since 
we  then  have  5  <  0.  As  seen  in  Fig.  3,  after  the  dip  in  the  radiation 
pattern  disappears  when  0  *  0O,  the  farfield  level  is  reduced  for  any  addi¬ 
tional  increase  in  <p. 

In  Eq .  (5),  it  is  the  denominator  of  the  directivity  function  which 
accounts  for  the  effect  of  variations  in  0  on  the  dip  in  the  radiation 
pattern.  As  a^  is  increased,  the  angular  dependence  of  the  denominator  is 
reduced,  making  the  aperture  factor  in  the  numerator  more  significant.  The 
aperture  factor  attains  its  maximum  value  of  unity  when  9V  *  0  and 
9X  *  (oji/co.) (J)  =  a,  where  a  is  the  direction  of  (see  Fig.  1).  Thus, 

for  sufficiently  high  attenuation  of  the  primaries,  Eq .  (5)  yields  results 
which  are  consistent  with  those  for  dispersive  plane  wave  interaction. ‘ 

From  Figs.  3  and  4,  we  find  that  the  radiation  pattern  is  not  shifted 
significantly  with  respect  to  the  z  axis  when  i>^0.  For  example,  when 
0  -  0O*  the  peak  is  shifted  less  than  one  degree  off  axis,  whereas  0*8.5°. 
This  is  because  we  assumed  low  attenuation  (aj*0.1)  to  emphasize  the 
effects  of  dispersion.  Such  low  attenuation  permits  the  interaction  region 
to  extend  beyond  the  near field  of  the  primaries,  whereby  conditions  do  not 
approximate  those  of  plane  wave  interaction. 
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